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In measurement-based quantum computation (MBQC), local adaptive measurements are per- 
formed on the quantum state of a lattice of qubits. Quantum gates are associated with a particular 
■ measurement sequence, and one way of viewing MBQC is that such a measurement sequence prepares 

f"^ ' a resource state suitable for 'gate teleportation'. We demonstrate how to quantify the performance 

(^i of quantum gates in MBQC by using correlation functions on the pre-measurement resource state. 
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A recent breakthrough in quantum computing has been the realization that quantum computation can proceed 
solely through single-qubit measurements on an appropriate quantum state [l], [3| ■ The canonical example of such a 
resource state is the cluster state [3, 0], which is a universal resource for MBQC on suitable lattices or graphs 0]. A 
handful of other universal resources for MBQC have recently been identified [3, HI, [H, 0] , but it is still not known what 
properties of a quantum many-body system allow for MBQC to proceed. 

In recent work Q, we have proposed that the ability to perform MBQC on a quantum many-body system may be 
identified using appropriate correlation functions as order parameters. This claim stems from the observation that 
MBQC is a means of preparing resource states for gate teleportation 0, 0, 03 ■ With a cluster state, it is possible by 
local measurements and feedforward alone to prepare such resource states allowing gate teleportation for a universal 
set of gates between essentially any set of qubits in the lattice. The performance of MBQC can be determined by 
calculating the fidelity between the resource state that is actually prepared and the ideal resource state. Here, we 
demonstrate how to express the resource state after the measurements in terms of correlation functions of the original 
J> ■ state prior to the measurements. These results provide an alternate perspective to Theorem 1 in Q, which shows 
0^ ! that the gates in the cluster state MBQC scheme function because of certain correlations in the original cluster state. 
■ In particular, our results apply to characterize gate performance in quantum states that are not the cluster state; 
with such correlation functions, it is possible to directly quantify the suitability of a given quantum many-body state 
^ . • for performing such MBQC gates. 

We note that our work has a close relation to the concept of localizable entanglement [ll[ . For the state of a quantum 
many-body system, the localizable entanglement between two arbitrary qubits is defined as the maximum amount 
of entanglement that can be created between these two qubits by performing local measurements on the remaining 
qubits. If this entangled state is viewed as a resource for quantum teleportation, then the localizable entanglement 
serves to quantify the ability to perform the trivial or 'identity' gate (i.e., teleportation) using local measurements. 
We note that the localizable entanglement in some systems can be quantified by correlation functions [ll|, [IH , using 
similar techniques as described here. Our work generalizes these results by considering non-trivial quantum gates, 
which include multi-qubit quantum gates. 

The paper is structured as follows. In Sec. [TT1 we review some of the essential terminology and mathematical 
structure of cluster-state quantum computing. Sec. IIIII presents the key general results of the paper, relating MBQC 
quantum gates to correlation functions. The correlation functions for the identity gate arc calculated explicitly in 
Sec. IIV1 and non-trivial gates including the 7r/2-gate, Hadamard, and Z-rotation gates are presented in Sec. El Two- 
dimensional gates are addressed in Sec. IVH and a general method for concatenating gates in Sec. IVIIl We finish with 
some brief conclusions in Sec. IVIII1 

II. PAULI OPERATORS, STABILIZERS, AND THE CLUSTER STATE 

The Pauli matrices are labeled X,Y,Z and are defined as 
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The group generated by the Pauli matrices under matrix multiplication is known as the Pauli group G±. The Pauli 
group G n on n-qubits is defined as G n = G% ® G\ ® • • • G\ . The Clifford group on n qubits is defined to be the 
group of unitary operators that map the Pauli group onto itself; i.e., U is in the Clifford group if UAU^ 1 <E G n for 
all AeG n . 



A. Stabilizer Formalism 



The stabilizer formalism is a method of describing a state of a quantum system by specifying a set of eigenvalue 
relations instead of its components in some basis. The eigenvalues of a complete set of commuting observables 
completely specifics a state. We define the stabilizers of a quantum state to be the set of operators for which the 
state is a +1 eigenstate. Clearly, any two stabilizers must commute, and thus the set of stabilizers forms an Abelian 
group. This group can be specified by its generators, and homomorphisms on the group can be specified completely 
by their effect on the generators. As a result it is sufficient to study the generators, a subset of the group, instead of 
the whole group. For example, the state 1 00) can be described as the state which is stabilized by the operators Z\ 
and Zi\ the stabilizer group of this state is generated by these two stabilizer operators. The stabilizer formalism was 
first used to describe quantum error correction codes [13], but are widely applicable to a variety of other situations. 

The standard stabilizer formalism is defined to only allow elements of the Pauli group G„ to be used as stabilizers. 
A set of Pauli stabilizers satisfies some key properties: 

• Elements of G n either commute or anti-commute. 

• If A, B, C £ G n and A anticommutes with both B and C, then A commutes with BC. 

The power of the stabilizer formalism lies in its ability to compactly describe certain quantum states, as well as 
their evolution under Clifford group operations and Pauli measurements. First, an arbitrary n-qubit state requires 
2 n complex numbers to completely describe, by specifying the contribution from each of the 2™ basis vectors of the 
n-qubit system. For systems which are stabilized by the Pauli group, a set of n stabilizers is sufficient to describe an 
n-qubit stabilizer state (l3| . i.e., n stabilizers can form a complete set of commuting observables. Hence if applicable, 
the stabilizer formalism offers a compact way of denoting states of a quantum system. Second, the stabilizer formalism 
is very efficient in describing a state under the evolution of unitary operators belonging to the Clifford group, the 
group of operators which map the Pauli group back to itself under conjugation. It is also efficient for describing 
the evolution of a state of a multi-qubit system under a projective measurement in the X, Y , or Z basis. A simple 
prescription exists which tells us how to obtain the stabilizers of the post-evolution/measurement state from the 



pre-evolution/measurement state 14 1. 



B. Cluster State 

The cluster state is a many-qubit entangled state Q • Consider an arbitrary graph with a qubit at each vertex; the 
cluster state on this graph is characterized by the stabilizers 

Ka = X a ® n£n bhd(a) Z n , (2) 

for each qubit a of the state, where nbhd(a) is the set of qubits adjacent to a (via the graph structure). In other 
words, the cluster state |$) satisfies K a \&) = |$) for all a. 

Two examples of cluster states on graphs will be considered here: a one-dimensional line, for which the stabilizers 
take the form Ki = Zi—iXiZi+i, and a two-dimensional square lattice, for which the stabilizers take the form 
Ki = XijZij-iZij+iZi-ijZi+ij. 



III. EXPRESSING GATES AS CORRELATION FUNCTIONS 



Consider a one-dimensional lattice of qubits prepared in the state po- Singling out two qubits, a- ln and a ou t, we 
wish to consider a measurement sequence on the remaining qubits in the lattice that yields a two-qubit resource state 
on qubits a; n and a ou t for gate telcportation. Let J label the measurement outcomes, and Pj be the corresponding 
projector. Following the measurements, a unitary Uj conditional on J is applied to qubit a ou t. Averaged over all 
possible measurement outcomes, the resulting resource state is 



P = EjUjPjPoPjU 



(3) 
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Equivalcntly, we can characterize this resource state using expectation values of bipartite Pauli operators A (g> B on 
qubits a in and a out , as 

(AB) = ZjT[[(AB)UjPjpoPjU$\ = E/&[(ABj)PjPoPj] , (4) 

where Bj = UjBUj. The set of such correlation functions, for A,B £ {I, X, Y, Z} spanning the set of Pauli operators, 
will completely specify the bi-partite resource state. 

For each of the MBQC gates given in 0, we can make use of a remarkable relation: there exists a string of operators 
S acting on some set of the measured qubits which is independent of the measurement outcomes, and an operator B' 
on a out that is also independent of J, such that 

(ABj)Pj = PjASB' . (5) 

Now, using the projector properties Pj = Pj and J2j Pj = ^ gives 

(AB) = Tr[(ASB')p ] . (6) 

Thus we can relate the resource state prepared after the sequence of measurements to a correlation function of 
the original state po prior to measurements. That is, the correlation functions characterize the pos^-measurement 
resource state using expectation values of strings of operators on the pre-measurement state. This argument is trivially 
extendible to multi-qubit gates. 

It is critical to this development that one can identify such a string of operators S. In the examples presented, 
the correction operator is a Pauli operator, and essentially it is the simple algebraic properties of the Pauli operators 
which arc responsible for the existence of S. That the corrections arc always Pauli operators makes the analysis of the 
Clifford gates especially simple, and we find that S for such operators takes the form of a product of Pauli operators. 
Even for our non-Clifford gate, we can still identify an appropriate operator (this time, a sum of product of Pauli 
operators). 



IV. THE IDENTITY GATE 



We will consider performing an identity gate between two qubits in a line with an odd number of qubits between, 
qubit k and qubit k + 21. To perform this gate one measures X on each of the 21 — 1 qubits between these two, 
and measures Z on qubit k — 1 and the qubit k + 21 + 1. We label measurement outcomes by the eigenvalues of the 
measured operators, which for Pauli operators are either +1 or — 1. Specifically, we label the measurement outcome 
for the X measurement on the qubit k + j by m,j = —1, 1, and the outcome for the Z measurement on the qubit k — 1 
(k + 2j + 1) by mi (m r ). We can then define two parities 



Pz = 



3=1-1 

1 - mim r ( ] [ m 2 j) 

3 = 1 



/2, 



Px = 



3=1 

i-(n™2,-i)]/2, 

3 = 1 



and write the correction unitary 



77 , — Y px 7 Vz 



(7) 



(8) 



We now seek to identify the commutation identity of Eq. ([5]). The unitaries only act on the B space (qubit k + 21) 
so operators that act only on the A space arc unchanged on commuting through the projectors and do not depend on 
the measurement outcome. For the B operators Xk+21 and Zk+21, we need to absorb the factors of ±1 that arise from 
the correction unitary by adding factors of the measured obscrvables. If we perform an X measurement for example 
then mP m = P m X, where P m is the projector on the m eigenstate. We have 



3=1-1 3 =l 

Xj = UjX k+ 2iU} = mim r ( ] [ m 2 j)X k+2 i , Zj = UjZ k+2 iU} = (JJ m 2j -i)Z k+2 i 

3=1 3=1 



(9) 



and therefore we have the relations 



3=1-1 3 =l 
XjPj = PjZ k ^x( ] [ Xk + 2j)Xk+2lZk+2l+l , ZjPj = Pj(Y\. Xk+2j-l)Zk+2l ■ 

3=1 3=1 



(10) 
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FIG. 1: Measurement sequences to create the resource state for gate teleportation of various single-qubit gates. The resource 
state qubits are aj n = 1 and a out — 5, and qubits 0, 2, 3, 4, and 6 are measured as shown. 



At this stage, we note that the correlation functions we would expect for an ideal identity gate are (XX) = (ZZ) 
-(YY) = 1. Using Eq. ©, we find 



3=1 

(XX) = Tr Z fc _i(JJX k+2j) Zk+2l+lP0 
3=0 
3=1 

(ZZ) = Tr Z k (Y\_ X k+2j -i)Z k+ 2i + ipo 

3 = 1 



3=1 

= Tr[(l[K k+2j )p 

3=0 
3=1 



Tr[(n#fc+2i-i)po 

3 = 1 



(11) 

(12) 



The correlation function (YY) follows similarly. As these correlations are equivalent to expectation values of cluster 
stabilizers, if po is the cluster state these expectation values will both be unity. 

Other expectation values, for example local expectation values (XI) or (IX), or those involving any other Pauli 
operator combination, can also be explicitly determined. The set of all such correlation functions will completely 
characterize the resource state. 



V. OTHER SINGLE-QUBIT GATES 



We now turn our attention to other single-qubit gates. Rather than directly consider correlation functions over 
arbitrary lengths, we restrict our attention to a fixed length (specifically, 3 intermediate qubits) between ai n and a out . 
In Section IVHi we will show how to concatenate such fixed-length gates together (possibly with the identity gate) to 
form gates of arbitrary lengths. Specifically, we consider creating bipartite resource states between qubits a- m = 1 and 
«out = 5 by measuring qubits 2, 3, and 4 along with Z measurements on qubits (to the left) and 6 (to the right). 
See Fig. [I] 



A. The Hadamard gate 



Here we construct the correlation functions for the Hadamard gate H = . We perform the Hadamard 

gate, by measuring Z on qubits and 6 and measuring Y on qubits 2, 3 and 4. This will give us correlation functions 
between qubits 1 and 5. The ideal resource state for Hadamard gate teleportation satisfies (XZ) = (ZX) = 1. 
Defining 



Pz 



/2. 



Px 



1 — mo 7713777,4 



/2, 



X^ x Z^ z 



1 — 777-2 TO3 777g 

the correction unitary for the Hadamard gate is then 

Uj-- 

We then have 

ZjPj = PjZoYzY^Zs 

and therefore the correlation functions are given by 

(XZ) = TrlZoXMYiZspo] = Tr[KxK 3 K m ] 

(ZX) = Ti[Z x Y 2 Y z X 6 Z tP o] = Tr[X 2 A' 3 ^ 5 Po] 

We note that these correlation functions can be expressed as expectation values of products of cluster stabilizers; they 
will yield a value of one if po is the cluster state. 



XjPj = PjY 2 Y 3 X 5 Z 6 , 



(13) 

(14) 
(15) 

(16) 
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B. The tt/2 gate 

The 7r/2 gate S = ^ is performed on a 1-D cluster state with qubit 1 as input and qubit 5 as output as in 

Fig. [U The gate is implemented by measuring Z on qubits and 6, Y on qubit 3 and X on qubits 2 and 4. The ideal 
resource state for 7r/2 gate teleportation satisfies (ZZ) = (X(—Y)) = 1. 
Defining 

Px 2 = t 1 - m 2m 4 ] /2 , p^ 2 = [l - mom 2 m 3 m 6 ] /2 , (f 7) 

the correction unitary for the n/2 gate is then 

Uj=xfx 2 zf /2 . (18) 

The relevant correlation functions are then 

(ZZ) = TcfaXiXtZtpo] = Tr[K 2 K iPo ] , 

(X(-Y)) = -Tv[Z X 1 Y 3 X 4 Y 5 Z 6 p Q ] = Tr [K^K^po] . (19) 

Again, we note that these correlation functions can be expressed as expectation values of products of cluster stabilizers; 
they will yield a value of one if po is the cluster state. 



C. A non-Clifford gate 

We now consider a non-Clifford gate - a rotation U z (9) by angle 8 about the Z axis. Again, we consider three 
intermediate qubits, as in Fig. [1] 

The ideal resource state p on A and B for gate teleportation of U z (6) satisfies (ZZ) = (XX-g) = 1, where 
X_g = U Z (—8)XU Z (9) = cos 9X — s'mdY. In MBQC, such a resource state is prepared on qubits 1 and 5 by 
measuring Z on qubits and 6, measuring X on qubits 2 and 4, and measuring X v = U z {rf)XU z {—rf) on qubit 3, 
where r) = m 2 8. We note that the measurement basis on qubit 3 depends explicitly on the outcome of the measurement 
on qubit 2. The correction unitary on qubit 3 is 

Uj = xf x zf z (20) 

with 

Px = [l - m 2 m i ] /2 , <p% = [l - m Q m 3 m 6 ] /2 . (21) 

For the X measurements on qubits 2 and 4, we can use m^Pj = PjX 2 and m^Pj = P7X4, and for qubits and 6 we 
can use moPj = PjZq and m^Pj = PjZq. For the X^ measurement yielding result 777,3, the situation is slightly more 
complicated, because rj = m 2 9. However, it is straightforward to show that 

m 3 Pj = P 7 (cos 6X 3 + sin 8X 2 Y 3 ) . (22) 

Thus, wc have 

ZjPj = PjX 2 X 4 Z 5 (23) 

XjPj = Pj(cos 6Z X 3 X 5 Z 6 + sin 9Z X 2 Y 3 X 5 Z 6 ) (24) 

YjPj = Pj(cos 9Z X 2 X 3 X 4 Y 5 Z 6 + sin eZ Y 3 X 4 Y 5 Z 6 ) . (25) 

Note that the right hand side of the equations is of the desired form: independent of the measurement results. 

These results allow us to express the two-qubit expectation values for the post-measurement resource state in terms 
of correlation functions on the pre-measuremcnt state. We have 

(ZZ) = T±\{Z x X 2 XaM = Tr[K 2 K iPo ] (26) 

(II-s) = Tr[(Z A 1 A 3 A 5 Z 6 (cos 2 8 + sin 2 8Z 3 X 4 Z 5 ) 

+ cos0sin0Z o A 1 A 2 y 3 A 5 Z 6 (l - Z 3 X 4 Z 5 ))p ] 

= Tr[(^iif 3 if 5 (cos 2 8 + sin 2 8K 4 ) + cos 6 sin 6(Z Q Y 1 Z 2 )K 2 K 3 {1 - K 4 )K 5 )p ] . (27) 

The term ZqY 4 Z 2 is not a cluster state stabilizer, and has an expectation value of for the cluster state. All other 
terms are stabilizers of the cluster state, and both correlation functions can be seen to have an expectation value 1 
on the perfect cluster state. In addition, this result agrees with Eq. (fT9"|) for 8 = tt/2. 
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VI. TWO-DIMENSIONAL CLUSTER STATES AND THE CSIGN GATE 



The single-qubit gate sequences, and their corresponding correlation functions, can be straightforwardly generalized 
to the cluster state on a two-dimensional square lattice. One-dimensional 'strips' can be created in the square lattice 
by performing Z measurements (and their corresponding Pauli corrections) to remove qubits from either side of the 
strip. 

However, for the purposes of defining simple correlation functions, it is easier to define single-qubit gate sequences 
along diagonal lines, as in Fig.[2ja). Such diagonals eliminate the need for Z measurements along the sides of the strip, 
and they are only required at the ends. Consider an example where we label qubits such that a m is at coordinate (f , 1) 
and a ou t is at coordinate (n,n). For an ideal cluster state, the products of stabilizers Jl'Z™ and JX-=i 1 ^i+M-i 
along and parallel to this diagonal arc themselves stabilizers. By measuring qubits (i, i) for 1 < i < n and qubits 
{i + l,i — 1) for 1 < i < n — l,wc use the standard rules for updating stabilizers to give 



fi-n 



i—n—l 



n * 



U=2 ""M 



i -r-ri=n— 1 

1 — 1L=i 



Z\ \Z n ,nZi aZ n +i 

,n— 1 



(28) 
(29) 



From this expression, we see that it is only necessary to measure qubits the six 'end' qubits in the Z-basis in order 
to obtain the two-qubit resource state stabilized by Xi^X n/n and Zi \Z n ^ n ; it is not necessary to measure the qubits 
on either side of the diagonal strip. 

With such diagonal strips, the correlation functions for the identity gate take the form 



j— n 

(XX) = Tr[([[K ili 



j—n — 1 

(ZZ)=Tt[{ [] Ki+M-i)A> 



(30) 



We consider the simplest version of a two-qubit gate: the CSIGN gate (a Clifford gate) that also implements a 
cross-over of control and target qubits. This gate is defined in terms of its action on Pauli operators as 



CSIGN 




L Oout 

Za mt ® X bo 

4ou t ® z bou 



(31) 



The measurement sequence is illustrated in Fig. [5J The relevant correlation functions from the above 'input-output' 
relations that will characterize the CSIGN gate are the expectations of four products of stabilizers: 



(X ain X aout Z bo J = Tr[K aiD K 3 K aoatPo ] , 
{Z ain X bin X bo J = Tr[K bin K 4 K bmtPo ] , 
(Z ain Z ao J=Tr[K 1 K 4 p \ , 
(Z bin Z bo J = Tr [K 2 K 3 po] , 



(32) 
(33) 
(34) 
(35) 



can be appended with diagonal strings of stabilizers in the direction of the arrows (and terminated with Z measure- 
ments as in Fig. Ufa)) to reach distant qubits. With X measurements on qubits 1-4, the resulting state provides the 
CSIGN transformation. 



VII. CONCATENATING GATES 



It is straightforward to concatenate many Clifford gates together into a single gate, and to calculate the resulting 
correlation functions. The essential idea is to equate the 'output' qubit of the first gate with the 'input' qubit of the 
second. For the gates sequences defined here, recall that Z measurements are performed at one qubit beyond each 
end of the gates. In this case where we combine two gates, the Z measurement prior to the 'input' end of the second 
gate and the 'output' qubit of the first gate are not performed, and their corrections are to be left out. Finally, an X 
measurement is performed on this joining qubit. 
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FIG. 2: (a) A measurement pattern on the cluster state that localizes entanglement between sites ai n and a ut, where X (Z) 
denotes a measurement in the X-basis (Z-basis). The two strings of stabilizer products, centred on sites connected by the 
parallel diagonal lines, directly quantify the fidelities of single-qubit gates between ai n and a ou t in MBQC. (b) The measurement 
sequence corresponding to the CSIGN gate between a and b. 



As we have expressed the correlation functions of our Clifford gates in terms of cluster stabilizers, it is straightforward 
to determine the correlation functions describing a combined gate: one simply takes the product of the corresponding 
stabilizer operators. 

We note that it is still possible, though less straightforward, to concatenate non-Clifford gates. The difficulty 
with non-Clifford gates is that they involve adaptive measurements, wherein the measurement basis can depend on 
the entire 'past history' of the computation. As a result, it is typically impossible to write the general form of the 
correlation functions for non-Clifford gates in a way that is independent of the specific choice of prior gates. However, 
in any given situation, such correlation functions can be defined using the methods presented here. 



VIII. CONCLUSIONS 



We have presented a general method for expressing the performance of quantum gates in the cluster-state model 
of MBQC as correlation functions on the pre-measurement resource state. With such correlation functions, viewed as 
order parameters, one can investigate the existence of a robust ordered phase in various models of quantum many-body 
systems that will allow for MBQC to occur. One such model is considered in [8[. 
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